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, ABSTRACT 

' We present the results of L-band spectroscopical observations of local bright Ultralu- 

I minous Infrared Galaxies (ULIRGs), performed with ISAAC at the VLT. The excellent 

, sensitivity of the telescope and of the instrument provided spectra of unprecedented 

■ quality for this class of objects, which allowed a detailed study of the AGN/starburst 

contribution to the energy output, and of the composition of the circumnuclear ab- 
ly^ . sorber. We discuss the L-band spectral features of seven single sources, and the statisti- 
' cal properties of a complete sample of 15 sources obtained combining our observations 
fH , with other published 3-4 ^m spectra. Our main results are: 

1. When a spectral indicator suggesting the presence of an AGN (low equivalent width 
of the 3.3 /im emission line, steep X — fx spectrum, presence of an absorption feature 
at 3.4 /im) is found, the AGN is always confirmed by independent analysis at other 
wavelengths. Conversely, in all known AGNs at least one of the above indicators is 



o 



^ \ present. 

2. Two new diagnostic diagrams are proposed combining the above indicators, in which 
starbursts and AGNs are clearly and completely separated. 
k>( , 3. The above diagnostic techniques are possible with spectra of relatively low quality, 

■ which can be obtained for several tens of ULIRGs with currently available telescopes. 

' This makes L-band spectroscopy the current best tool to disentangle AGNs and star- 

bursts contributions in ULIRGs. 

4. The L-band properties of ULIRGs are heterogeneous. However, we show that all 
the spectral differences among ULIRGs can be reproduced starting from pure intrinsic 
AGN an starburst spectra and varying two parameters: the amount of dust extinction 
of the AGN component, and the relative AGN/starburst contribution to the bolomet- 
ric luminosity. 

5. Using the above decomposition model, we show that AGNs in ULIRGs have a low 
dust-to-gas ratio and a dust extinction curve different from Galactic. 

6. The estimate of the presence and contribution of AGN in a complete sample show 
that AGN are hosted by '--^ 2/3 of ULIRGs, but their energetic contribution is relevant 
(> 30% of the bolometric luminosity) only in ~ 20% of the sample. 

Key viTords: galaxies: active - galaxies: starburst ~ infrared: galaxies 



1 INTRODUCTION 

Ultraluminous Infrared Galaxies (ULIRGs, Lm > 10^^ Lq) 
are the most luminous sources in the local Universe. More- 

* Based on observations collected at the European Southern Ob- 
servatory, Chile (proposals ESO 69.A-0643, ESO 73.B-0574). 



over, they represent the local counterparts of the class of 
high redshift objects dominating the far-infrared (FIR) and 
sub-mm backgrounds (Sanders & Mirabel|4^ and references 
therein). Understanding the origin of the huge infrared emis- 
sion of ULIRGs is therefore important in itself, given the rel- 
evance of this class of objects in the low redshift Universe, 
and in order to understand the nature of the high redshift 
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FIR sources, which are too faint to be studied in detail with 
current instrumentation (even if the Spitzer Space Telescope 
is a major improvement in this direction, Werner et al. 15711 . 

These considerations prompted a large interest in ob- 
servational and theoretical studies on this class of sources. 
The main well established aspects on the physics of ULIRGs 
are: (a) the infrared (from a few /xm to ~ 1000 /^m) emission 
dominates the total luminosity, and is due to dust reprocess- 
ing of higher frequency radiation (b) ULIRGs are mainly 
found in interacting or highly irregular systems (Sanders et 
al.lii). 

Two different physical processes are possible as the pri- 
mary energy source of ULIRGs: nuclear activity due to gas 
accretion onto a supermassive black hole and/or strong star- 
burst activity. Several studies have been performed in the 
past 20 years, with the aim of understanding (a) which frac- 
tion of ULIRGs host an AGN, and (b) which is the contri- 
bution of the AGN to the bolometric luminosity. 

In the optical and near-IR, spectral classification of 
ULIRGs of the IRAS 1 Jy sample (Kim & Sanders 13) es- 
tablished that only a minor fraction of ULIRGs have a clear 
signature of an AGN (Veilleux et al. i5i Veilleux, Kim, & 
Sanders Is^. Veilleux, Sanders, & Kim l54) . This fraction in- 
creases with luminosity from ~ 20% for objects with lu- 
minosity in the range Lir ~ 10^^ — 10^^ Lq to ~ 50% 
for objects with Lir > 10^^'^ Lq (Veilleux et al. Is^) . 
A large fraction (~ 40%) of ULIRGs in the 1 Jy sample 
have optical/near-IR spectra typical of Low Ionization Nar- 
row Emission Regions (LINERs). This classification leaves 
the question about the starburst/AGN relative contribution 
open. 

In the mid-IR, ISO spectra of bright ULIRGs provided 
a powerful diagnostics for the study of the nature of these 
sources (Genzel et al. [ll hereafter G98). The presence of 
high ionization coronal lines in the mid-IR has been proven 
to be an indicator of the presence of an AGN (G98, Lutz 
et al. Isol) . while the presence of strong emission features 
due to polycyclic aromatic hydrocarbon (PAH) molecules is 
an indication of starburst dominance (see next Section for 
further details). 

In the X-rays, spectra of bright ULIRGs, obtained with 
XMM-Newton by our group (Franceschini et al. fl^ Braito 
et al. 0, Braito et al. |^ and with Chandra (Ptak et al. Issl) 
revealed the capability of hard X-ray spectroscopy to disen- 
tangle the contributions of starbursts and AGNs. 

The works briefly referred to above signiflcantly en- 
hanced our understanding of the nature of the emission of 
ULIRGs, but still have relevant limitations, which can be 
divided into two classes: 

(a) Reliability of the diagnostic techniques. Most of the AGN 
indicators fail in several cases to detect an existing AGN, if 
this is heavily obscured by gas and dust. An extreme exam- 
ple of such limitations is the ULIRG NGC 6240, which has 
no indication of the presence of an AGN at any wavelength 
with the exception of the hard X-rays, where the AGN direct 
emission escapes the high column density {Nh > 10^'* cm~^) 



^ In this paper we will not discuss powerful radio loud quasars, 
which have an infrared luminosity in the ULIRG range, mainly 
due to synchrotron emission of relativistic particles in the jet, but 
are not dominated by infrared radiation. 



absorber at i5 > 10 keV (Vignati et al. l55i) . and the flat re- 
flected emission, with a prominent iron Ka emission line 
with equivalent width EW > 1 keV are observed in the 2- 
10 keV range. 

(b) Extension to faint sources. The most powerful indica- 
tors of the presence of buried AGNs in ULIRGs have been 
proven to be the hard X-ray spectra, and ISO mid-IR spec- 
tra. In both cases a good signal-to-noise spectrum is needed 
in order to detect the AGN/starburst indicators: in the X- 
rays, a high S/N is required to disentangle the hard power 
law due to the AGN emission from the thermal component 
due to the starburst; in the mid-IR a high S/N is needed to 
have an accurate determination of the continuum and, con- 
sequently, a reliable measurement of the equivalent width of 
the PAH emission features. With current infrared detectors, 
the required spectral quality can be achieved only for a small 
number (~ 10 — 15) of bright objects. A big improvement 
for the mid-IR analysis is expected thanks to the Spitzer 
Space Telescope, which will be able to provide high quality 
mid-IR spectra for several tens of ULIRGS, as first results 
have already confirmed (Armus et al.^. 

An interesting new way of studying ULIRGs has been 
proven to be L-band (~ 3—4 fim) spectroscopy. Imanishi and 
collaborators (in particular, Imanishi & Dudley 2000) dis- 
cussed the physics of the emission/absorption features in the 
3-4 ^m interval, and showed that low-resolution spectra of 
the brightest ULIRGs, obtained at 4 meter class telescopes, 
provide effective diagnostic elements in order to disentangle 
the AGN and starburst contributions to the observed emis- 
sion. Recently, Imanishi, Dudley & Maloney (2005) applied 
this method to investigate the energy source in a sample of 
~ 40 objects from the 1 Jy sample (Kim & Sanders 1998). 

Risaliti et al. (2003) presented a high signal-to-noise L- 
band spectrum of the bright ULIRG IRAS 19254-7245 ob- 
tained with the instrument ISAAC at the Very Large Tele- 
scope (VLT) in Paranal, Chile. This high quality spectrum 
suggested that 8 meter class telescopes are able to provide 
good enough spectra for (a) a detailed study of the phys- 
ical properties of bright ULIRGs, and (b) an analysis of 
relatively faint ULIRGs, for which other indicators are not 
available with the current instrumentation. 

In this paper we present new L-band spectra of bright 
ULIRGs obtained with the instrument ISAAC at the VLT 
in Paranal, Chile, and we perform a complete analysis of 
the L-band diagnostics applied to a representative sample 
of ULIRGs, consisting of our new observations and previ- 
ously available L-band spectra. In Section 2 we review the 
L-band spectral features useful in order to investigate the 
energy source in ULIRGs. In Section 3 we present our sam- 
ple. In Section 4 we present the new L-band spectra. The 
physical analysis of each new spectrum is presented in Sec- 
tion 5. In Section 6 we present a detailed discussion of our 
results, and an analysis of the different indicators in a com- 
plete flux-limited sample, obtained merging our sample of 
southern sources with other spectra obtained by Imanishi 
and collaborators for northern ULIRGs. In particular: (1) 
we present new two-dimensional diagnostic diagrams where 
starburst and AGNs are completely separated (2) we discuss 
the physical properties of the dusty circumnuclear medium 
in ULIRGs hosting AGNs, showing that a different amount 
of hot dust absorption can explain the main differences 
among spectra of these sources; (3) we show that the amount 
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of hot dust can also explain the differences among starburst- 
dominated ULIRGs. 

Our conclusions are summarized in Section 7. Through 
this paper we estimate luminosity distances using the 
WMAP concordance cosmology {ho — 70, fim = 0.3, fix = 
0.7, e.g. Spergel et al.lsol 



2 L-BAND DIAGNOSTICS FOR ULIRGS 

Disentangling the contribution of heavily obscured AGNs 
and starbursts to the luminosity of ULIRGs is intrinsically 
complicated, because of the reprocessing, in both cases, of 
the primary optical/UV radiation into mid- and far-IR emis- 
sion. 

However, several physical differences between the two 
energy sources can have direct consequences on the observed 
reprocessed radiation in the infrared: 

• The X-ray emission of an AGN, much stronger than 
that of a starburst (with respect to the bolometric emis- 
sion) affects the gas and dust composition around the central 
source. For example, (a) highly ionized atomic species, re- 
quiring a strong ionizing X-ray continuum, are present only 
around AGNs, and (b) complex hydrocarbon dust grains can 
be present only around starbursts, while they cannot survive 
in the close vicinity of AGNs, since they are destroyed by 
the X-ray photons. 

• An active nucleus is an extremely compact radiation 
source with respect to a starburst region. This implies that 
strong absorption effects in AGN spectra are easily obtained 
through compact and dense absorbers along the line of sight. 
The same effects are harder to obtain in starbursts, since the 
absorber should be large enough to cover the whole emitting 
region. 

• The amount of high temperature (from a few 10"^ K 
to the sublimation temperature, Tsub '~ 1, 500 K) dust is 
higher in AGNs than in starbursts. This affects the infrared 
spectral energy distribution, in particular in the near- and 
mid-IR. 

In particular, the emission in the L-band is due to re- 
processing of the primary radiation by hot dust (temper- 
ature in the range 100-1000 K), i.e. by material relatively 
close to the central source. As a consequence, the observed 
L-band spectra are heavily affected by the above physical 
differences. 

Here we briefly review the main diagnostic features in 
ULIRGs, most of which are present in some of the spectra 
shown in Fig. 1. These indicators have been calibrated on 
objects for which the starburst/ AGN classification is known 
from other, independent studies. However, several cases are 
known for which some of these diagnostic tools fail. We 
will discuss these cases, and possible improvements of the 
diagnostic techniques, in Section 6. 

3.3/im emission feature. A strong, broad emission 
feature at a rest-frame wavelength A — 3.3/im is the most 
prominent non-continuum feature in the spectra of most 
ULIRGs. The origin of this emission feature is believed 
to be the reprocessing of UV radiation by polycyclic aro- 
matic hydrocarbon (PAH) molecules. These complex organic 



molecules (for a review on their structure and physical prop- 
erties, see AUamandola, Tielens, & Barker^ are responsible 
of several broad emission features in the mid-IR wavelength 
band, and have been extensively studied in ISO spectra of 
ULIRGs (e.g. G98). 

The 3.3 /xm emission feature is particularly interesting 
in this context because it is an isolated feature, which can 
be easily disentangled from the continuum emission. This is 
not the case in many mid-IR ISO spectra, where the combi- 
nation of broad emission features due to PAHs with broad 
absorption features due to dust makes the determination of 
the continuum level (and therefore of the equivalent width 
of the lines) extremely uncertain (one notable exception is 
the PAH emission feature at 6.2 jj,m, as discussed in Fischer 
©)• 

Both theoretical and observational evidence suggests 
that the strength of the 3.3/im PAH feature is an indicator 
of the relative contribution of AGN and starburst to the 
emission of ULIRGs: 

- Observationally, it is found that objects which are inde- 
pendently classified as starburst-dominated show a 3.3 
feature with an equivalent width EW3.3 '^lOO nm or higher. 
This is true both at ULIRGs luminosities and for objects 
with luminosities of the order of ^ 10^^ Lq (Moorwood l34l 
Imanishi & Dudlev ITsl) . Objects known to be dominated 
by AGN emission show a weak or absent {EW3.3 < 30 nm) 
3.3 fim emission feature (Imanishi & Dudley ,15,). As we will 
discuss in more detail later, this is not a "perfect" indicator. 
However, the correlation is strong, and no exception is 
known at least in cases of extreme values (i.e. all known 
sources with EWs.i < 50 nm or EWi.3 > 100 nm are 
dominated by AGN or starburst, respectively). 

- From a physical point of view, two mechanisms are 
expected to be effective in decreasing the equivalent width 
of the 3.3pm line in AGNs: (a) the X-ray radiation can 
destroy the PAH molecules which are not shielded by 
a column density from a few lO'^^ to 10^* cm~^ of gas, 

tending on the geometrical and physical conditions (Voit 
, and (b) the strong L-band continuum (much stronger 
than in starbursts due to the presence of a larger amount 
of hot dust) can dilute the 3.3/im emission from outer or 
gas-shielded starburst regions. 

absorption feature. A deep absorption fea- 
ture is clearly present in the spectrum of several ULIRGs. 
The most prominent cases are IRAS 19254-7245S (Fig. 1, 
see Risaliti et al. |43for details), UGC 5101 (Imanishi et al. 
113) and IRAS 08572-^3915 (Imanishi & Dudley [S), where 
the optical depth^ at 3.4 ^m is T3.4 ~ 0.8 — 1. The origin 
of this feature is believed to be absorption by hydrocarbon 
dust grains, in particular by vibrational stretching levels of 
the C-H2 and C-H3 groups (Sandford et al. and refer- 
ences therein, Pendleton & AUamandola [s^)- In the high- 
est signal-to-noise detections, as the ones mentioned above, 
a detailed study of the chemistry of the absorbing dust is 



^ All the optical depth estimates presented in this paper arc com- 
puted as the natural logarithmic ratio between the estimated con- 
tinuum flux and the measured emission at the wavelength corre- 
sponding to the mimimum of the absorption profile. 
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possible through the analysis of the substructure_of the ab- 
sorption profile (Risaliti et al. l42l Mason et al. .'i.'l i. 

The presence of the 3.4/im absorption feature is an 
indicator of a strong AGN contribution to the ULIRG emis- 
sion. Again, the motivation has both a physical explanation 
and an observational confirmation: 

- Observationally, all the above mentioned ULIRGs are 
independently known to host a powerful AGN. The other 
known cases among ULIRGs where this feature has been 
revealed (Imanishi & Dudley even if with a smaller 
optical depth, are all classified as AGNs on the basis of in- 
dependent classification at other wavelengths. Furthermore, 
L-band spectra of non-ULIRG obscured AGNs show strong 
3.4/xm absorption flmanishi ild) . 

- Physically, in order to detect a high 3.4/im absorption 
optical depth, a high column density absorber in front of 
the primary source is required. Imanishi & Maloney (19) 
showed that dilution and saturation effects prevent the 
observed optical depth from being higher than ~ 0.2 in 
case the absorbing dust is spatially mixed with the energy 
sources. Therefore, in order to achieve high values of t^a, 
an external screen completely covering the radiation source 
is needed. Since the absorbing dust is expected to lie in gas 
clouds, we can estimate the amount of gas expected to be 
associated with the dust assuming a Galactic dust-to-gas 
ratio, and the extinction curve of Pendleton et al. (1361) . 
We find JVh ~ 3 - 5 x 10^^T3.4 cm"^. The actual value 
could easily be more than an order of magnitude greater, 
given that the measured dust-to-gas ratio in AGNs is 
systematically lower than Galactic (Alonso-Herrero, Ward, 
& Kotilainen 0, Maccacaro, Perola, & Elvis Maiolino et 
al. l3lh . Such high column densities are easily present around 
AGNs (~half AGNs in the local Universe are absorbed by 
column densities higher than 10^'* cm~'^, Risaliti, Maiolino, 
& Salvati ,4(1') , while a complete covering of a starburst 
region would imply a huge amount of gas: assuming the 
linear dimension of a starburst region emitting ~ lO^'^ Lq 
to be of the order of or greater than one kiloparsec, the re- 
quired mass of gas would be Mgas > 10^^ x (_R/100 pc)^ PIq. 

S.l^m absorption feature. A broad absorption fea- 
ture at ~ 3.1 has been detected in several ULIRGs 
(Imanishi & Maloney ITgl) . The origin of this feature is 
absorption by ice-covered dust grains. Analogously to the 
3.4/im feature, it can be shown that optical depths rs.i > 0.3 
are not possible unless an external screen covering the source 
is present, if a fraction of 30% of dust grains covered by ice is 
assumed (as estimated for the nearby starburst M82, Iman- 
ishi & Maloney 2003). Then the same argument as above 
suggests that this is highly unlikely in starburst-dominated 
sources. This expectation is confirmed by the discovery of 
broad ice absorption in AGN-dominated ULIRGs (Imanishi 
& Malonev Il9ft . However, the uncertainty on the possible 
profiles and widths of this absorption features makes it less 
straightforward than the previous indicators, especially in 
low signal-to-noise spectra. In particular, in many cases it 
is not possible to distinguish in AGN-dominated objects 
whether a concave spectrum in the 3-3.6/im region is due 
to a broad ice absorption feature, or to a change in spectral 
slope (in a A — /a plane) between the near-IR region (still 
dominated by the tail of the direct emission of the AGN 
accretion disc) and the mid-IR region (dominated by repro- 



cessing of hot dust). In order to distinguish between these 
scenarios, a high signal-to-noise spectrum is required, in or- 
der to constrain the continuum slope using the feature-free 
3.5-4 wavelength region. Alternatively a broader wave- 
length coverage, involving K-band and/or M-band spectra 
are needed. Moreover, a different (higher) fraction of dust 
grains covered by ice could increase the upper limit of rs.i 
in starbursts. 

In the next Section we will separately discuss each case 
among our spectra where ice absorption could be present. 

Continuum slope. The observed continuum slope in 
the L-band spectra of ULIRGs has a wide range, spanning 
from r =~ —2 to F 3 when fitted with a fx — 
model. Observations of known AGNs and starbursts suggest 
that positive F are an indication of AGN activity. The 
physical explanation of this trend is clear: a high F implies 
strong dust reddening of a compact source, and hence the 
presence of an AGN. As for the other indicators, we note 
that observations confirm the effectiveness of the continuum 
slope in extreme cases (F > 2 are always associated to 
AGNs) , while objects with F < 2 (i. e. most of the ULIRGs 
observed so far) can be either AGNs or starbursts. We will 
go into further details on this point in Section 6. 

Atomic emission lines. Atomic emission lines are in 
principle powerful diagnostics of the presence of an AGN. A 
broad atomic emission line, such as the Bra line at 4.05 fim, 
or a high ionization line such as the [SilX] A3. 95 fim or the 
[MgVIII] A3. 028 fj,m, are strong indicators of the presence 
of an AGN. These diagnostics have been successfully used 
on nearby Seyfert 2 Galaxies (Lutz et al. lisl) . 

However, out of these three lines, by far the most in- 
tense expected in the (rest frame) L-band, two {Bra and 
[SilX]) are in all our objects out of the well calibrated 2.9- 
4.1 fj.m interval, due to the redshift effect, and the third 
one, [MgVIII] is in the poor atmospheric transmission inter- 
val (Fig. 1). Therefore, we will not use this diagnostics in 
this work. 



3 THE BRIGHT ULIRGS SAMPLE 

The sample chosen for L-band spectroscopic observations 
is the one of G98, consisting of 15 ULIRGs in the IRAS 
Bright Galaxy Sample (BGS, consisting of all the galaxies in 
the IRAS all sky survey with galactic latitude ||6|| > 5° and 
60 flux density 5*60 > 5.24 Jy, Sanders et al. 46), with two 
variations: one source (IRAS 23060-1-0305) has been dropped 
because it does not fulfill the flux requirement of the BGS, 
and two sources (IRAS 08572+3915 and IRAS 05189-2524) 
have been added, because they have 5*60 > 5.24 Jy (Table 1). 

The luminosity criterion used to select ULIRGs is the 
standard one, i.e. Lm > 10^^ Lq, where Lm is the total 
8-lOOO^m luminosity, calculated using the IR flux equation 
of (|4a): 

Fir = 1.8 X 10-''(13.485i2 + 5.125*25 + 2.58560 + 5ioo) (1) 

where _Fi_r is the total IR flux in units of erg s~^ cm~^, and 
5i2, 525, 560, 5ioo are the flux densities in the four IRAS 
filters, in units of Jy. 



L-band spectroscopy of ULIRGs 5 



One source, NGC 6240, has a slightly lower total in- 
frared luminosity. However, it is included in our sample, as 
well as in that of G98, because all the morphological and 
physical properties appear to be typical of ULIRGs. 

Two different considerations drove us in choosing this 
sample: 

1. These are the IR brightest ULIRGs, therefore we expect 
to obtain the best possible quality for our L-band spectra. 
This is crucial at the present level of development of L-band 
diagnostics: as briefly discussed in the introduction, we still 
need to perform detailed studies of bright objects in order 
to better understand and calibrate the physical diagnostics, 
before applying this technique to a larger sample of fainter 
sources. 

2. The sources in our sample are the best studied ULIRGs 
at other wavelengths. In particular, complete near-IR pho- 
tometry and spectroscopy and ISO mid-IR spectroscopy 
is available, and in the X-rays we recently performed a 
complete X-ray spectroscopic study with XMM-Newtqn 
(Franceschini et al. Braito et al. Braito et al. Q). 
This will allow us (a) to investigate multi-wavelength indica- 
tors of the AGN/starburst contribution to the IR emission, 
and (b) to compare our results in the L-band with exist- 
ing AGN/starburst classifications based on observations at 
other wavelengths. 

Out of the 16 ULIRGs in our ULIRGs sample, we 
selected eight sources visible from Chile, and with no 
previous L-band spectroscopic observations. For six more 
sources (IRAS 05189-2524, MKN 231, MKN 273, UGC 5101, 
Arp 220, IRAS 08572-1-3915) in the sample 3-4/xm spectra 
obtained at the 3.6m UKIRT and 3.0m IRTF telescopes have 
been published by Imanishi & Dudley Jisl) and Imanishi, 
Dudley, & Maloney (H^ . We will include their results in the 
discussion in Sect. 4. 

For one source (NGC 6240) with an L-band spectrum 
available in the literature (Imanishi & Dudley HH) , we re- 
cently obtained a new, higher quality spectrum with ISAAC. 
The results of this observations will be published elsewhere 
(Risaliti et al. 2005, in prep.). In the following we will refer 
to the results of this new observation. 

For one source (IRAS 19254-7245S) we already pub- 
lished a detailed spectral analysis (Risaliti et al. 142") . The 
results of this analysis will be included for completeness in 
tables and plots, but will not be further discussed. 

We report in Tab. 1 the basic known photometric data 
of our sources. 



4 OBSERVATIONS, AND DATA REDUCTION 

The L-band spectra were obtained with the instrument 
ISAAC at the UTl unit of the Very Large Telescope (VLT) 
on Cerro Paranal, Chile. All but two objects were observed 
during three nights on June 02-04, 2002. IRAS 20100-4156 
and IRAS 22491-1808 were observed on July 30, 2004. 

The sky conditions were photometric, with a seeing of 
the order of 1 arcsec in K (0.7 arcsec for IRAS 20100-4156). 
For all but one observations a 20x1 arcsec slit was used. 
For one source (IRAS 20100-4156) we used the 20x0.7 slit. 
We note that the seeing in the L band is typically better 
than in K. This is confirmed by the profile analysis made 
on our observations (see below for details) , which show that 



in all cases the fraction of slit loss is lower than 10%. 6 
out of 8 sources (the exceptions being IRAS 20551-4250 and 
IRAS 22491-1808) are known mergers, with two distinct nu- 
clei resolved in infrared imaging. For all these objects we 
chose the slit angle in order to obtain spectra for both nuclei. 
The angles were chosen based on 2MASS images, and on the 
short images taken before each spectral observation. In three 
cases (IRAS 12112+0305, IRAS 14348-1447, IRAS 23128- 
5919) we were able to obtain high quality spectra of both 
nuclei. In one object (IRAS 19254-7245) we obtained a high 
quality spectrum for the brightest nucleus, and a low S/N 
spectrum for the faintest one, enough for an estimate of the 
flux and continuum slope. In two objects (IRAS 17208-0014 
and IRAS 20100-4156) we obtained a spectrum of the bright- 
est nucleus only. Out of the two single nucleus sources, we 
obtained a high quality spectrum of IRAS 20551-4250, while 
IRAS 22491-1808 was not detected. 

The on-source time varies from one hour for the bright- 
est objects to ~ 3 hours for the faintest. The observations 
were performed with the low resolution grating (LR) in the 
L band, with a coverage of the ~ 2.9 — 4.1 /im wavelength 
band with a resolution A/AA ~ 360. 

The main data regarding the observations are listed in 
Table 2. 

In order to avoid saturation due to the high background 
(~ 3.9 mag per arcsec^ in L-band) all the spectra were 
taken in chopping mode, with single exposures of 0.56 s. 
The spectra were then aligned and merged into a single im- 
age. We performed a standard data reduction, consisting 
of flat-fielding, background subtraction, and spectra extrac- 
tion, using the IRAF 2.11 package. 

For each object we also acquired the spectrum of a stan- 
dard star immediately before or after the target observation. 
The star was chosen in order to have the same airmass as 
the target, within 0.1, and with an L magnitude between 4 
and 6. This allowed us to obtain high quality spectra with 
2 minute long observations. 

Corrections for sky absorption and instrumental re- 
sponse were obtained from the standard stars spectra. All 
the standards are B stars, therefore we assumed a pure 
Raleigh- Jeans intrinsic spectrum in the L-band. 

In order to facilitate the background subtraction, the 
observations were performed using a mosaic consisting in a 
shift of the slit by 20 arcsec up and down along the chosen 
direction. However, the instrumental response is not con- 
stant along the slit direction (this effect is mostly removed 
by the fiat-fielding; however, given the high level of noise, 
residual differences would strongly affect the final result). 
We compensated for this effect performing a separate re- 
duction for each of the three spectra relative to the three 
different positions along the slit. The results before the di- 
visions by the three standard star spectra (obtained in the 
same way) were in general quite different one from the other, 
both in absolute fiux and in spectral shape. However, after 
the division by the stars, the three spectra of all sources 
overlapped within the errors. This gives us an important 
self-consistency check for our spectral reduction procedure. 

In order to obtain a precise absolute calibration, we 
took into account aperture effects by analyzing the profiles 
of both the targets and the standard stars along the slit. We 
assumed a Gaussian profile and we estimated the fraction of 
fiux inside the slit assuming a perfect centering. We estimate 
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Name 


z 






oa 


ca 
^100 


^IR 








orf 
^5.9 


■^7.7 


TRAS 1 21 1 


0.072 


0.06 


0.29 


5 


5.8 


1.1 






14.3 


5.9 


7.7 


TR A 1 91 1 9-Ln'^nf^N"P 
iJri,/\o IZllZ-t-UoUOiN-Cj 






n 91 


o.o 


A 9 


U.O 












IRAS 14348-1447S 


0.082 


0.07 


0.33 


4.5 


4.7 


1.3 


15.6 


14.7 


13.7 


1.3 


1.9 


IRAS 14348-1447N 




0.03 


0.17 


2.4 


2.4 


0.7 


16.1 


15.1 


14.3 






IRAS 17208-0014A 


0.043 


0.2 


1.7 


31.1 


34.9 


2.3 


14.0 


12.9 


12.4 


3.2 


4.9 


IRAS 19254-7245S 


0.062 


0.19 


1.15 


5.3 


5.6 


1.1 


14.0 


12.7 


11.7 


7.0 


9.5 


IRAS 19254-7245N 




0.01 


0.05 


0.2 


0.2 


0.2 












IRAS 20100-4156 


0.130 


0.1 


0.3 


5.2 


5.2 


3.6 


15.0 


14.3 


13.7 


1.4 


2.1 


IRAS 20551-4250 


0.043 


0.3 


1.9 


12.8 


10.0 


1.1 


13.4 


12.2 


12.2 


4.5 


7.0 


IRAS 22491-1808 


0.078 


0.1 


0.6 


5.4 


4.5 


1.5 


14.7 


14.0 


13.6 


1.5 


4.4 


IRAS 23128-5919N 


0.044 


0.05 


0.3 


2.9 


2.7 


0.2 


14.7 


13.7 


13.3 


4.2 


6.4 


IRAS 23128-5919S 


71 


0.15 


1.3 


7.9 


7.3 


7.0 


14.3 


13.4 


12.6 







Table 1. Photometric NIR to FIR data for our sources. In five cases we report NIR photometry for two components of the sources. 
These are the cases where we were able to disentangle the two components in our L-band observations. For these objects we scale the total 
IRAS fluxes according to the K magnitudes of the single nuclei. Notes: °': flux density in the IRAS filters, in units of Jy; Total infrared 
luminosity in units of 10^^ Lq, estimated from Eq. 1 and assuming the concordance cosmology (H(,,Q,m,^\) = (70,0.3,0.7), Spergel 
et al. l50l) . Infrared magnitudes from D ue, M irabel, & Maza J^, except for the K magnitudes of the two nuclei o f IR AS 12112+0305, 
which are from Kim, Veilleux, & Sanders i2j[) ; Mid-IR flux densities obtained with ISO, from Rigopoulou et al. i39l) . 



Source Observation Date T"^ Slit width Pos. Angle*" Standard Star Type L'^ 



IRAS 12112-1-0305 


2002-06-02 


180 


1" 


35 


Hip 076126 


B3V 


5.92 


IRAS 14348-1447 


2002-06-03 


180 


1" 


25 


Hip 065630 


B3IV 


6.52 


IRAS 17208-0014 


2002-06-02 


60 


1" 





Hip 076126 


B3V 


5.92 


IRAS 19254-7245/1 


2002-06-02 


60 


1" 


170 


Hip 076126 


B3V 


5.92 


IRAS 19254-7245/2 


2002-06-03 


60 


1" 


170 


Hip 000183 


B4V 


5.42 


IRAS 20100-4156 


2004-07-30 


90 


0.7" 





Hip 011407 


B5IV 


4.65 


IRAS 20551-4250 


2002-06-02 


60 


1" 





Hip 076126 


B3V 


5.92 


IRAS 22491-1808 


2002-06-02 


60 


1" 





Hip 076126 


B3V 


5.92 


IRAS 23128-5919/1 


2002-06-02 


60 


1" 





Hip 038593 


B2V 


5.98 


IRAS 23128-5919/2 


2002-06-03 


60 


1" 





Hip 000183 


B4V 


5.42 



Table 2. Observation log of our ISAAC program. 2 sources were observed twice. Notes: ": observing time in minutes. *": Slit position 
angle, with respect to the north-south inclination. The angle increases in counter-clockwise direction, starting from the north. '^: L 
magnitude of the standard star. 



this procedure to have an error of the order of ^ 10 — 15% 
or better. In three cases (Table 2) this is confirmed by the 
comparison of the spectra obtained in two different nights, 
and corrected with two different standard stars: the final 
cahbrated spectra are always consistent within the errors. 
A further consistency check is provided by the comparison 
with the magnitudes at longer wavelengths (when available 
from ISO observations), which show a good agreement with 
our results. 

For each source the final spectrum, obtained (a) merg- 
ing the three single spectra relative to the three slit posi- 
tions, and (b), when two or more observations were present, 
merging the final calibrated spectra, was rebinned in order 
to have a statistical significance of each bin of at least 2a. 

The error bars were estimated from the Poissonian noise 
in the sky counts, which are by far the dominant source of 
noise (for ease of comparison, our brightest source has a 
magnitude L~ll, i.e. a ~ 600 times lower fiux than the sky 
background in the same extraction region). 

The final L-band spectra of our sources are shown in 
Fig. 1. 



5 DATA ANALYSIS 

As a first step for the interpretation of the L-band spec- 
tra of our sample of ULIRGs, we fitted the "feature-free" 
continuum with a single power law. The chosen wavelength 
interval is the 3.4-4 fim (rest frame) for all sources ex- 
cept IRAS 19254-7245S and IRAS 20551-4250, where we 
took into account the conspicuous absorption features at 
~ 3.4 fj,m. The best fit power laws are plotted in Fig. 1, in 
order to better analyze the absorption and emission features. 

As a second step, we fitted the whole spectra in order to 
quantitatively measure the relevant parameter of the major 
spectral features. The best fit model consists of a power law, 
an emission line at rest frame wavelength Xrest ~ 3.3 fim, 
and when needed an absorption feature (with Gaussian pro- 
file) at Xrest ~ 3.4 fim (hydrocarbon dust grain absorption) 
and one at Xrest ~ 3.1 fim (water ice absorption). The re- 
sults of the spectral fits are shown in Tab. 3. 

Here we briefiy discuss the spectrum of each source, 
estimating the AGN/starburst contribution and comparing 
our results with those obtained through the observations at 
other wavelengths. In particular we refer to G98 for ISO 
mid-IR spectroscopy, and to Franceschini et al. ilSl) for 
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Figure 1. Flux-calibrated L-band spectra of our sample of ULIRGs. Wavelengths are in the rest frame. The shaded vertical band shows 
the spectral region with bad atmospheric transmission. In each spectrum, the continuous line represents a power law model fitted on the 
spectral region red-ward of the broad emission feature at 3.3/im (i.e. at ~ A > 3.4/im in the rest frame). Exceptions are IRAS 19254- 7245S 
and IRAS 20551-4250, where the two obvious absorption features at ^ 3.4 — 3.5/im were also excluded from the fitting region. 



XMM-Newton X-ray spectroscopy. 

IRAS 12112+0305. The two nuclei in this source 
have similar L-band spectra. In both cases the high equiv- 
alent width of the 3.3 fim line, the absence of 3.4 
absorption and the fiat spectrum suggest a starburst origin 



for their infrared luminosity. This is in agreement with the 
classification at other wavelengths: 

- in the optical, the [OIII]/If/3 ratio, and the low flux of 
[Nil] A6583A, [SII] A6716A and [OI] A6300A with respect 
to Ha, suggest a LINER classification (Veilleux, Kim, & 
Sanders Isa). 
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Figure 1. - Continued. 



- in the mid-IR the high equivalent width of the 7.7 fim PAH 
emission feature suggests a pure starburst classification. 

- in the hard X-rays no indication of AGN activity is found. 

IRAS 14348-1447. The L-band spectrum of both 
nuclei is dominated by a high equivalent width 3.3 /im 
emission feature (EW>100 nm). The continuum is ~ flat 
in both cases, with an uncertainty in the brightest nucleus 
(S) due to a possible presence of a moderate {t ^ 0.1 — 0.2) 



ice absorption. If we allow for ice absorption, the best fit 
continuum slope is F ~ —0.2, while without ice absorption 
r ~ 0.5. The values in Tab. 3 refer to the fit to a power 
law+emission line model, without absorption features. This 
model provides a good representation of the continuum, 
with a small excess in the 2.9-3.1 /im, which could be 
due to the tail of the K-band emission coming from a 
broader region than the L-band emission. The starburst 
classification is confirmed by ISO mid-IR spectroscopy and 
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name 


pa 


3.3 


tS , 
3.1 


3.4 




oLV Ayo.o/ -^(^J 


IRAS 12112+0305SW 


-0.64±0.3 


78±21 


0.5±0.1 


- 


13.4 


9.5 




-i.4itU.O 


y i ±04 


U.OitU.l 




io.o 


y.o 


IRAS 14348-1447S 


0.55±0.3 


98±23 






12.7 


9.9 


IRAS 14348-1447N 


-0.13±0.6 


100±37 






13.5 


9.6 


IRAS 17208-0014 


-0.32±0.1 


83±5 


0.4±0.05 




11.7 


9.7 


IRAS 19254-7245S 


2.30±0.03 


26±3 




0.8±0.05 


10.3 


10.6 


IRAS 19254-7245N 


0.1±1.2 








13.9 


9.2 


IRAS 20100-4156 


0.85±0.7 


92±41 






13.3 


10.0 


IRAS 20551-4250 


2.75±0.13 


72±9 




0.3±0.1 


12.0 


9.6 


IRAS 23128-5919S 


-0.2±0.1 


43±5 


0.45±0.05 




11.1 


10.0 


IRAS 23128-5919N 


-1.05±0.2 


130±22 






12.7 


9.4 



Table 3. Spectral analysis results. Slope of the continuum power law, defined as oc A'". '': Equivalent width of the 3.3 /^m emission 
feature, in units of nm. Optical depth of the 3.1 /^m ice absorption feature. Optical depth of the 3.4 /xm hydrocarbon absorption 
feature. L magnitude. ^ : luminosity Xlx at 3.5 /^m, in log of bolometric solar luminosity. 



X-ray spectroscopy. 

IRAS 17208-0014. The L-band spectrum is dom- 
inated by a strong PAH 3.3 /im emission feature. The 
continuum slope F ~ 0.3 is indicative of a starburst origin 
of the L-band emission. An ice absorption feature at 3.1 /jm 
is clearly detected. A single Gaussian absorption component 
is significant at a 5 cr statistical level. This, together with 
IRAS 23128-5919S, is the best constrained ice absorption 
feature among the sources in our sample, thanks to the high 
S /N spectrum which allows a precise continuum determina- 
tion. The estimated optical depth is r3,i = 0.4 ±0.05. In the 
mid-IR and X-rays IRAS 17208-0014 appears dominated 
by starburst activity. 

IRAS 19254-7245. The high quality spectrum of the 
southern nucleus show both a 3.3 /im emission feature with 
an equivalent width (EW'~25 nm), and two clear AGN 
features, i.e. a steep continuum (F = 2.3) and a strong 
absorption at 3.4 ^m (t = 0.8). The substructures of the 
absorption profile can be used to constrain the chemical 
composition of the absorbing dust. A detailed analysis 
of this spectrum is presented in Risaliti et al. Jl^)- The 
northern nucleus is extremely faint, and a detailed spectral 
analysis is not possible. The continuum is fiat, and no 
spectral feature is present. The low S/N prevents us from 
performing any further analysis. A detailed study of the 
multiwavelength spectral energy distribution led to the 
conclusion that the infrared emission of this nucleus is due 
to an inactive galactic component 10). 

IRAS 20100-4156. Only the brightest nucleus of this 
source is detected. The spectrum (the lowest S/N among 
the primary nuclei in our sources) is dominated by the 
3.3 /im emission line, suggesting a starburst origin of the 
L-band emission. The continuum slope is quite uncertain, 
due to the low S/N. An argument similar to that for 
IRAS 12112-0305NE can be applied: the continuum can 
be fitted (a) by a flat power law, with a residual excess 
at 2.9-3 /xm, or (b) by a steeper (F ~ —0.5) continuum 
plus an ice absorption feature with r3.i ~ 0.5. In both 
interpretations, no clear indication of the presence of an 
AGN is found, even if we cannot rule out a possible minor 
AGN contribution, due to the low statistics available. 



The mid-IR spectrum of this source suggests a starburst 
classification. In the X-rays, no clear indication of an AGN 
is found, however the location of this source in the X-ray 
diagnostic diagrams of Franceschini et al. is between 
the typical starbursts and AGNs positions, suggesting a 
possible contribution of an AGN. 

IRAS 20551-4250. The high quality spectrum of this 
source is the most puzzling in our sample. A strong 3.3 yLxn 
emission feature (EW=70 nm) suggests a dominance of the 
starburst contribution. However, a steep and inverted con- 
tinuum slope (F ~ 2.7) and the presence of a 3.4 /xm absorp- 
tion feature (r3.4 ~ 0.3) points to the presence of an AGN 
dominating the continuum emission. This is a clear exam- 
ple where the simple classification, based on the spectral 
diagnostics described in the previous section, fails. We will 
discuss this source further in Section 6, where an improve- 
ment of the above indicators will be discussed. The mid-IR 
spectrum of IRAS 20551-4250 shows clear indications of the 
presence of an AGN. Similarly, the X-spectrum is typical of 
a heavily absorbed AGN-dominated source. The presence of 
the AGN was also revealed by the analysis of the polarized 
optical spectrum Jj^l . 

Risaliti et al. i42l) and Mason et al. showed that a 
detailed investigation of the chemistry of the dust grains re- 
sponsible of the 3.4 iim absorption is possible, based on the 
analysis of the substructure of the 3.4 fj,m absorption feature. 
Such analysis on the only other source of our sample show- 
ing the 3.4 ^m absorption feature, IRAS 19254-7245, led 
to the conclusion that the dust composition in that source 
is remarkably similar to that observed in Galactic stellar 
sources. 

The continuum-subtracted spectrum of IRAS 20551- 
4250 in the 3.3-3.8 yLxn range is plotted in Fig. 2. A 
hydrocarbon dust absorption feature is clearly visible in 
Fig. 2, and confirmed by the fit with a single Gaussian 
absorption component, which is significant at > 3 cr. 
A detailed analysis of the chemical composition of the 
absorbing dust (analogous to that done in R03 and Mason 
et al. (33) is not possible here, because of the low S/N and 
the low optical depth (r ~ 0.3 when fitted with a single 
Gaussian profile). However, from a qualitative analysis of 
Fig. 2 we can tentatively note a stronger absorption dip 
at ~ 3.5 ^m with respect to the Galactic source IRS7 
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Figure 2. Continuum-subtracted spectrum of IRAS 20551-4250 
in the 3.3-3.8 ixra wavelength range, com pared with the spectrum 
of the Galactic center source IRS 7 i36l) . The complex structure 
of the hydrocarbon dust absorption feature is clearly visible. The 
first dip is due to vibrational transitions of pure aliphatic hydro- 
carbon molecules; the second dip is due to the same transitions 
in hydrocarbon molecules where electronegative group, such as 
-NO2 or -OH, are present. 



(Pendleton et al. |3fl|) and IRAS 19254-7245 (Risaliti et 
al. li^l . This would imply that the hydrocarbon molecules 
responsible for absorption in this spectral region are richer 
than in the other two sources of electronegative groups such 
as -NO2 and -OH. 

IRAS 23128-1808. The southern nucleus has a flat 
continuum, with a moderate equivalent width 3.3 /xm 
emission (EW3.3 = 43 ± 5, well below the typical value for 
starbursts). A strong ice absorption (T3.1 = 0.45) is clearly 
detected. The relatively low EW of the 3.3 ^m suggests the 
presence of an AGN. The 3.1 /xm absorption feature also 
suggests the presence of a compact source. In this case, 
the good continuum determination leaves no doubts on the 
interpretation of the 2.9-3.3 spectrum, since the alter- 
native scenario of an inverted continuum and a 2.9-3 pm 
excess is ruled out by the data. The northern nucleus has 
a typical starburst spectrum with a steep continuum and 
a strong 3.3 /xm emission line (EW~ 130 nm). No AGN 
indication is found in mid-IR spectra, while the hard X-ray 
spectrum clearly reveals the presence of an AGN. This 
source is therefore a case where the L-band analysis proves 
to be more effective than ISO spectroscopy in revealing the 
presence of an active nucleus as the origin of the infrared 
emission. 



Summary. We analyzed the L-band spectra of 7 
ULIRGs, 6 of which consisting of two distinct nuclei. All 
the single nuclei were detected except the faintest one in 
IRAS 20100-4156. The main results can be summarized as 
follows: 

• A direct evidence for the presence of an AGN (equiva- 
lent width of the 3.3 /xm emission feature EW3.3 < 50 nm, 
and/or steep and inverted continuum with F < — 1, and/or 
3.4 yLxn absorption feature) has been detected in three cases 
(IRAS 19254-7245, IRAS 20551-4250, IRAS 23128-5919). 
All these three sources have been classified as AGNs in our 
XMM-Newton X-ray survey of ULIRGs (Franceschini et al. 
d. 

• In two of the above cases the AGN indicators are not 
in agreement: the continuum of IRAS 23128-5919S is not 
inverted, and EW3.3 of IRAS 20551-4250 is typical of a pure 
starburst. 

• In two more nuclei (IRAS 12112-f 0305SW and 
IRAS 17208-0014) the detection of a strong ice absorption 
feature at 3.1 /xm is unambiguous. No other AGN indica- 
tion is present in these two sources either in the L-band or 
at other wavelengths. 

• A possible absorption feature in the 3.0-3.2 ^m range 
is detected in four more nuclei (IRAS 12112-f 0305NE, 
IRAS 14348-1447S and N, IRAS 20100-4156). In these cases, 
the interpretation as ice 3.1 iim absorption is not unique. 
An alternative fit can be obtained with a fiatter or inverted 
continuum, plus a 2.9-3.0 /im excess due to the tail of the 
near-IR emission of the host galaxy. 



6 ANALYSIS OF A COMPLETE SAMPLE 

The L-band spectra described above represent a significant 
step forward in our study of ULIRGs, for three main reasons: 

• High quality. The high signal-to-noise allowed a pre- 
cise determination of the parameters of the main spectral 
features in most spectra (Tab. 3). 

• Completeness. Our small sample consists of the 
sources of the G98 sample visible from Chile at the time 
of the observations. Since the G98 sample is made by the 
brightest known ULIRGs, and we did not introduce any 
physical bias, our sources can be considered as representa- 
tive of ULIRGs in the local Universe. Moreover, combi ning 
our results with those of Imanishi & Dudley (|fl3) and liTv 
we have good quality L-band spectroscopy of 13 objects out 
of the 15 of the G98 sample. This almost complete sample 
can be used to obtain statistical estimates on the average 
properties of ULIRGs in the local Universe. 

• Broad spectral coverage. The availability of obser- 
vations at other wavelengths (optical, near-IR, mid-IR, X- 
rays) allows a direct comparison of our results with those 
obtained with other, independent spectral diagnostics. This 
is fundamental in two respects: (a) it allows a precise cali- 
bration of the L-band starburst/ AGN indicators, and (b) a 
study of the complete spectral energy distribution of these 
sources is possible. 

We will use these properties in order to discuss three 
fundamental issues about the nature of ULIRGs: (1) Are 
L-band indicators capable to disentangle the AGN and star- 
burst contributions in all cases ? And is L-band spectroscopy 
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Name EWl;, log[(Ai^)3.5/i©]'' logLg.s/i^ L_ 



UGC 5101 


5+0.5 


35+10 


0.65+0.05 


10.1 


8.1 


0.95 


Arp 220 


0.05+0.05 


82+10 




9.3 


7.6 


1.5 


NGC 6240 


-0.1+0.1 


50+10 




9.9 
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10.9 
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0.014 
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-0.2+0.2 
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0.12+0.01 


11.3 


< 7.6 


0.58 


IRAS 20460+1925 


-0.3+0.2 


< 0.5 




11.7 


< 7.9 


3.3 


IRAS 23060+0505 


0.1+0.2 


< 1 




11.8 


< 8.3 


2.8 



Tab le 4. Spectral parameters of the ULIRGs in the com ple te sample and of five cont rol sources (in the last five lines). Data are from 
except for UGC 5101 (Imanishi, Dudley, & Maloney |3) , MKN 463 (Imanishi [H) , and NGC 6240 (Risaliti et al.2005). Spectral 
index of a power law fitting the 3.4-4.0 fivn continuum. The values are inferred from a visual inspection of the published spectra. 
Equivalent width of the 3.3 /im emission featur e. T he errors are conservative estimates based on visual inspection, except for the source 
UGC 5101, for which the error is reported in JlTt) . optical depth of the 3.4 fim absorption feature, "^i Log of the Xfx luminosity at 
3.5 fim, in units of Lq. Logarithmic ratio between the luminosity of the 3.3 /^m emission feature and the total IR luminosity. ^ : Total 
infrared luminosity in units of 10^^ Lq, estimated from IRAS fluxes according to Eq. 1. 



effective in finding AGNs in ULIRGs, with respect to indi- 
cators at other wavelengths? (2) How common are AGNs 
among ULIRGs ? (3) When an AGN is present, is it pos- 
sible to estimate the relative contribution of the AGN and 
starburst components to the total luminosity ? 

In the following, when the sample completeness is 
required (as, for example, in estimating the fraction of 
ULIRGs hosting an AGN) we will use the complete sam- 
ple described above. When completeness is not required, we 
will add to the above sample a few more sources, with lumi- 
nosity in the ULIRG range, or slightly lower, well known to 
be either starburst- or AGN-dominated in the L-band. This 
will increase the statistics and provide safer references for 
our interpretation. The relevant properties of these "control 
sources" are listed in Table 4. 

6.1 L-band Diagnostic Diagrams 

All the sources in our sample are well studied at other wave- 
lengths, and the starburst/ AGN contribution to their lumi- 
nosity has been deeply investigated. The availability of an 
independent classification is useful to calibrate our L-band 
indicators, and to check their effectiveness. In Tab. 5 we re- 
port the indications about the presence of an AGN from 
several L-band indicators, and from studies at other wave- 
lengths. 

The most important result emerging from the analysis 
of Tab. 3-5 is that in all cases where one indicator suggests 
the presence of an AGN (i.e. low EW3.3, T3.4 > 0, or F > 1) 
this is confirmed by independent classification. 

As a further step in the search for empirical indicators 
of the AGN/starburst presence in ULIRGs, we combined 
the main AGN/SB indicators in order to study the possible 
correlations among them. 

a) We plotted our sources in a EW3.3 — T plane 
(Fig. 3a). This plot shows that an almost 100% effective 
L-band classification is possible: AGNs and starbursts lie in 
two different regions in the EW3.3 — T plane, 
b) Fig. 3b shows that the same clear division is obtained 



by plotting EWs.s versus the R'f^^^ , defined as the loga- 
rithmic L-band to total infrared ratio, corrected for L-band 
continuum absorption. In order to estimate the absorption 
in the L-band, we made use of the relation Al ~ 12t3.4 
(Pendleton et al.^. 

These two results show the power of L-band spec- 
troscopy in detecting AGNs among ULIRGs. 

It is worth noting that the importance of the above 
analysis is related to the possibility of a large extension of 
the present work. 

If we only limit our study to the sample analyzed here, 
our results are only a confirmation of the conclusions already 
obtained through studies at other wavelengths, in particular 
X-rays. However, the relatively low signal to noise of the X- 
ray spectra of these objects (which are the brightest known 
ULIRGs) obtained with Chandra (Ptak et al.^ and XMM- 
Newton (Franceschini et al. 113) clearly indicates that it is 
not possible to extend the X-ray analysis to fainter objects. 
On the contrary, the quality of our spectra show that it is 
possible to obtain a rough estimate of the continuum slope 
and the 3.3 /im PAH feature equivalent width for a source 
as faint as L ~ 14. 

Considering the average L-band / FIR ratio (Fig. 3) , at 
least half of the sources in the IRAS 1 Jy sample of ULIRGs 
(Kim & Sanders 13) ^-re expected to be brighter than this 
limit. This implies that our method, that has been calibrated 
here thanks to the availability of multi-wavelength observa- 
tions of our bright sources, is capable to provide a statis- 
tical analysis on the presence of AGNs in a representative 
sample of 70 objects with redshift z < 0.15. (a first sam- 
ple of ~ 40 sources observed with Subaru (Imanishi, Dudley 
& Maloney 2005) is already available for this purpose, and 
more sources are scheduled for observation with ISAAC at 
the VLT in the next few months. This will be the subject of 
a future paper). 

Even more optimistic predictions are possible if we take 
into account that adaptive optics instruments are now avail- 
able for L-band spectroscopy, like the NAOS-CONICA at 
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Table 5. Detection of an AGN (either dominant or weak compared with the starburst component) in the complete sam ple of ULIRGs, 
according to several L-band, optical, mid-IR and X- ray spectroscopy. Re fere nces: 1: Imanishi , D udley, & M alone v Il7t 2: Imanishi & 
Dudley list 3; This work; 4: Veilleux, Kim, fc Sanders IsM 5: Veilleux et al.|53 6: Sanders et al.@ 7: Due, Mirabel, & Maza|3 8: Lutz 
Veilleux, & Genzel |28t 9: Imanishi et al.ES 10: Iwasawa et al. 12,4 11: Vignati et al. l55t 12: Prancescbini et al. Il.^t 13: Risaliti et al. l43t 
14: Laureijs et al. l26t 15: Severgnini et al. 143 
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Figure 3. Left: F — EW3.3 plot for our sample of ULIRGs and control sources. Right: equivalent width of the 3.3 fim emission feature 
versus the logarithmic ratio R'j°'~^ between absorption corrected L-band continuum luminosity and total IR luminosity. Errors on R'j°^^ 
are of the order of, or smaller than, the size of the points. The straight lines are empirical divisions between AGNs and starbursts. 



VLT. Since the noise in L-band observations of our sources 
is completely dominated by the thermal background (~ 
3.9 mag/arcsec^) a factor 10 smaller slit aperture (0.1 arc- 
sec instead of 1 arcsec) would imply a significant increase 
of the signal-to-noise, even when efficiency losses due to the 
adaptive system are taken into account. 

6.2 Physical properties of the nuclei of ULIRGs 

In the previous Section we discussed the effectiveness of the 
L-band indicators in detecting AGNs in ULIRGs. Even if we 
demonstrated that such indicators are indeed effective from 
a practical point of view, it is not obvious to find a relation 



between the variety of the observed L-band properties of our 
sample of ULIRGs and their physical properties. 

Here we want to investigate this fundamental issue, 
through a physical interpretation of the indicators discussed 
above and of the results of observations at other wave- 
lengths. 

Among AGNs, an apparent spread in spectral proper- 
ties suggests a high heterogeneity in the physical properties 
of the central emitting regions. However, two further 
indications suggest a simple physical interpretation: 
1) The correlation between the continuum slope, F, and 
the optical depth at 3.4 ^m, T3.4, is shown in Fig. 4. The 
sources with F > 1 are those with a high 3.4 fim absorption 
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Figure 4. 3.4 /^m absorption versus continuum slope for our sam- 
ple of ULIRGs and control sources. Typical errors on r3.4 are of 
the order of 10%. 



burst /AGN emission in ULIRGs will be discussed in the 
next Section 6.3. 

With regard to starburst dominated sources, the ob- 
served spread in the L-band spectral parameters (continuum 
slope r, EW3.3) is much smaller than in sources known to 
host an AGN. This spread could be either intrinsic (due to 
the age of the stellar populations) or due to an AGN com- 
ponent which has not been revealed through our indicators. 
For example, an unobscured or moderately obscured AGN 
contributing ~ 1/3 of the L-band emission would probably be 
missed by our analysis, and would imply a change by ~ 30% 
of the measurements of the above spectral parameters. We 
note that in order to have such a case the contribution of 
the AGN to the bolometric luminosity has to be negligible 
(of the order of a percent or smaller). 

Finally, we cannot exclude that a powerful AGN, even 
dominating the bolometric luminosity, with a high enough 
obscuration to make negligible its contribution in the L- 
band, is present in some of the starburst-classified ULIRGs. 
In this case the only affected parameter would be the L-band 
to bolometric ratio. Given the large spread in this parame- 
ter, (Fig. 4) even among starburst-classified sources alone, 
this possibility cannot be excluded for any of the sources in 
our sample. 



(t3.4 > 0.3). All the other sources are clustered around the 
(r ~ 0, EWx:i ~ 0) point in Fig. 3a. 

2) AU AGNs in the (F ~ 0, EW3.3 ~ 0) group (Fig. 3) are 
not heavily absorbed in the X-rays (column densities in the 
range 22 < logiVn < 23.5), while all the high F, high T3.4 
objects have log Nh > 24. Assuming a Galactic dust to 
gas ratio, log Nh ~ 23 implies tIl ~ 3. If the dust to gas 
ratio is lower than Galactic, as suggested in all cases where 
a direct measurement is possible (Maiolino et al. Isif) . the 
absorption in the L-band is even smaller. 

Combining the above indications, a physical interpreta- 
tion clearly emerges: 

- (F ~ 0, EW3.3 ~ 0) AGNs (Fig. 3) are not heavily ob- 
scured in the L-band. As a consequence, no strong absorp- 
tion feature is present, and the continuum emission is due 
to the direct emission of the hot circumnuclear dust heated 
by the optical/UV primary AGN emission. In these cases, 
the AGN component dominates in the L band even when a 
significant starburst contribution to the bolometric luminos- 
ity is present. This is due to the ~two order of magnitudes 
higher L-band/bolometric ratio in AGN than in starbursts. 
As a consequence, the PAH emission feature at 3.3 /im is 
either weak or absent. 

- (F > 1,T3.4 > 0.3) AGNs are heavily obscured in the L- 
band. The observed continuum is the result of a strong dust 
extinction by dust. The heavy dust obscuration results into 
a highly reddened continuum (high F) and a strong absorp- 
tion feature at 3.4 /^m. Since the AGN continuum is heavily 
absorbed, if a starburst is contributing to the bolometric 
luminosity, it is possible that it also represents a significant 
fraction of the L-band emission, depending on the amount of 
AGN extinction. This explains the wide range of measured 
values of the equivalent width of the 3.3 fim emission feature 
in these objects. 

A more quantitative treatment of the relative star- 



6.3 The AGN contribution to ULIRGs 

Estimating the relative contribution of AGNs to the lumi- 
nosity of ULIRGs is a more ambitious goal than the sim- 
ple detection of AGNs inside ULIRGs. A multi- wavelength 
analysis, based on models of the whole optical to FIR Spec- 
tral Energy Distribution (SED) is probably the most pow- 
erful approach to this problem. Such a study for several of 
our sources is currently being finalized (Fritz et al. 2005, in 
prep.). 

Here we only focus on L-band diagnostics. 

Our analysis of the AGN contribution to ULIRGs lu- 
minosity is done in two steps: (a) estimate of the frac- 
tion of ULIRGs hosting an AGN; (b) estimate of the rel- 
ative AGN/starburst contribution to the total luminosity in 
sources hosting an AGN. 

6.3.1 Fraction of ULIRGs hosting AGNs 

Out of the 14 ULIRGs in Table 5, only 5 do not show 
any evidence for the presence of an AGN. These sources 
belong to the flux-limited sample defined by the criterion 
Seo > 5.24 Jy (Sect. 2). In order to complete this sam- 
ple only two objects are missing: IRAS 15250-1-3609 and 
IRAS 22491-1808. The first one has never been observed 
in the L-band, while the second one was not detected in 
our ISAAC observation. Both sources have been observed 
in the X-rays with XMM-Newton (Franceschini et al. 
and in the mid-IR with ISO (G98), and no indication of 
the presence of an active nucleus has been found in both 
cases. Therefore, in order to be conservative in the estimate 
of the fraction of AGN-hosting ULIRGs, we classify these 
two sources as starburst-dominated. The final sample of 16 
ULIRGs consists of 7 starburst dominated objects (44%) 
and 9 objects with a significant AGN contribution (56%). 

These results have been obtained by using a flux selec- 
tion at 60 /im. This is not an unbiased selection criterion. 
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because of the higher fraction of bolometric luminosity emit- 
ted at 60 /xm by starbursts with respect to AGNs (see also 
the discussion in Risaliti et al. l4 J) . In our sample the average 
fraction of the total IR luminosity emitted at 60 ^tm (using 
Eq. 1) is 0.62 for the seven starbursts, and 0.52 for the 9 
AGNs. This implies that in order to have the same total IR 
flux limit for AGNs and starbursts, the flux limit at 60 
must be lower for AGNs by a factor (0.62/0.52=1.19). As- 
suming a euclidean logN-logS, the expected correction for 
the number of AGNs in a 60 /im flux limited sample is 
1.19^'''^ = 1.30. When this correction is taken into account, 
the final estimate of the fraction of ULIRGs hosting an AGN 
is 63%. 



6.3.2 Deconvolution of the AGN and starburst 
components 

An accurate estimate of the relative starburst/ AGN con- 
tribution for each source in the L-band requires a spectral 
decomposition of the two components. This is extremely dif- 
ficult due to the spread of intrinsic properties shown in Fig. 3 
and 4. 

A zeroth-order indicator of the relative contribution 
of the two components is the logarithmic ratio R'j^^^ be- 
tween the absorption-corrected continuum luminosity in the 
L-band and the total IR luminosity (Fig. 3b) . It is shown in 
Fig. 3b that sources whose L-band spectrum is dominated 
by the unobscured emission of an AGN have > —1, 

while the average value for starburst-dominated sources is 
j^corr ^ _2.6. For comparison , R l ~ —0.57 in the quasar 
spectral energy distribution of Hw . The huge difference be- 
tween the two values of R'^^^^ for starbursts and AGNs, 
and the relatively large dispersion around the average val- 
ues makes an accurate estimate of the relative contributions 
of the two components extremely difficult. For example, if 
50% of the total IR luminosity of a ULIRG hosting an un- 
obscured AGN is due to a circumnuclear starburst, R'^^^ 
would be changed only by a factor of two, well within the 
spread shown in Fig. 3b. The conclusion is that the L-band 
spectrum of ULIRGs hosting an unobscured AGN is dom- 
inated by the AGN emission even if the starburst is the 
dominant energy source. 

The situation is significantly different for obscured 
AGNs. In these cases dust extinction of a few magnitudes 
can make the AGN emission comparable to that of a star- 
burst of similar bolometric luminosity. The contiimum ex- 
tinction can be estimated from the measured 3.4-4 /im slope. 

In order to quantitatively estimate these contributions, 
we made several assumptions on the intrinsic AGN and star- 
burst emission, and developed a simple model to reproduce 
the observed L-band spectra. The main features of the model 
are the following: 

• The observed L-band spectra are the composition of 
AGN and starburst contributions, which are parametrized 
as: 



fx ~ afx^AGNC 



+ (1 - a)fx^sB 



(2) 



where a is the fraction of the L-band luminosity due to the 
AGN, and tl(A) is the wavelength dependent continuum op- 
tical depth in the L-band. fx^AGN and fx^sB are the intrinsic 
AGN and starburst L-band spectra. 



• Based on the results of the spectral analysis shown 
in Table 3 and 4, and in Figure 3 and 4, we assume 
fx,AGN = A""-^ and fx,SB = A"" '^ -I- fpAH. fpAH is the 
contribution of the 3.3 /xm emission feature, normalized in 
order to have EW-j,.z,sb ~ 110 nm with respect to the 
pure starburst continuum. The choice of the spectral index 
(r = —0.5 for both pure AGN and pure starburst) has been 
driven by the average values observed in our sample, both 
for starburst-dominated and AGN-dominated ULIRGs. The 
value for AGNs is in agreement with L-band spectra of pure 
type 1 AGNs (Imanishi & Wadalilh. The value for star- 
bursts is the average from our data, and has a large scatter 
(Fig. 3a). However, this uncertainty is not a major problem 
for our model (see details below). We note that our pure 
starburst spectra differ on average from the average spectra 
of lower luminosity normal galaxies (which have F ~ — 2, Lu 
et al. 27) , probably due to a larger contribution by hot dust 
in powerful starbursts. 

• We assume an extinction law tl{\) cx A~^'^^ (0). 

The model has two free parameters: tl and a, where tl ~ 
tl{X = 3.5 ^m). For each pair of these parameters we com- 
pute the slope, F, and the PAH feature equivalent width, 
EW3.3, for the composite spectrum. A good analytic ap- 
proximation (within a few percent) of these relations can 
be obtained noting that the wavelength interval (Ai,A2) 
used for the continuum fit is rather narrow (Ai ~ 3.4 /im, 
A2 ~ 4 ^m) . Using a first order expansion in 5 — (A2 — Ai)/Ai 
we easily obtain: 



EW3.: 



{l-a)EWx3.SB 



(1 



-0.5 + 



1.75Tae-^^ -I- (1 - a) 
ae-^L -f (1 _ Q,) 



(3) 



(4) 



We plot the results in Fig. 5. The two groups of curves 
are obtained freezing one of the two parameters and vary- 
ing the other. Therefore, they form a grid which can be 
directly used to estimate the relative AGN contribution and 
the AGN L-band extinction of each single ULIRG. 

Equations 3 and 4 can be easily solved in terms of a 
and Tl. A quantitative analysis of the results is shown in 
the Appendix. Here we concentrate on a visual analysis of 
Fig. 5, from which both the strengths and the limitations of 
the model clearly emerge: 

- It is possible to estimate with good accuracy the AGN 
contribution to the L-band emission in a large region of the 
parameter space (a > 0.1, rz, < 4). 

- Our diagnostics fails around the pure starburst region, 
where all the model lines converge and high degeneracy 
among the parameters is present. Therefore, we cannot ex- 
clude that L-band starburst-dominated objects host a heav- 
ily buried {tl > 5) AGN. However, in the cases where the 
AGN dominates, it is still possible to disentangle its contri- 
bution even in highly obscured objects, such as IRAS 20551- 
4250 and UGC 5101 (numbers 7 and 8 in Fig. 5). 

- The dependence on the choice of the intrinsic AGN and SB 
spectra is not strong: if, for example, the "pure starburst" 
point is moved (provided that it remains inside the starburst 
region of the plot), the only affected region is the one close 
to the point itself, where the model lines converge and the 
estimate of the parameters is not reliable. These limitations 
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Figure 5. Estimates of the relative AGN/SB contribution in a L-band continuum slope versus EW3.3 plot, according to the model 
described in the text. Dashed black lines represent combinations of AGN and SB with a fixed relative contribution with a varying L-band 
absorption of the AGN component (from = to tj^ = 6). The two numbers below each line are the fractions (in percent) of the AGN 
contribution to the L-band (first number) and to the bolometric emission (second number). Continuous red lines represent combinations 
with a fixed AGN absorption (as labeled for each line) and a varying relative contribution. Following a dashed line from the bottom-left, 
it starts from a no-absorption point and moves to the upper-right direction as absorption (i.e. reddening AND attenuation of the AGN 
component) increases. It then moves to the lower right direction because as the reddening further increases, the starburst contribution 
becomes more important than the absorbed AGN component. Following a continuous line, it starts from a right point representing a 
pure AGN, and then moves to the bottom-left as the starburst contribution increases. 



Numbers identifying each source refer to Table 1 



are quantitatively discussed in the Appendix, through the 
analysis of the statistical errors in the model. 

6.3.3 AGN/SB contributions to the bolometric luminosity 

A further step can be made estimating the AGN contribu- 
tion to the bolometric luminosity. This is possible taking into 
account the difTerent L-band to bolometric corrections for 
AGN and starbursts, Rl{AGN) and Rl{SB). We introduce 
a new parameter, K, defined as the ratio between the two 
bolometric corrections From the considerations made above 
{Rl{SB) ~ -2.6,Rl{AGN) ~ -0.6) we have K ~ 100. 

The fraction of the bolometric luminosity due to the 
AGN, aBOL, is then given by 

^BOL = -T (5) 

Q + A (i — Qj 



We list the results of the (a,asoL,T_L) estimates in 
Table 6. Restricting our analysis on the complete sample 
(i.e. neglecting the five "control sources") and adding up 
the contributions of both nuclei in the double sources, we 
obtain that the relative AGN contribution to the bolometric 
luminosity, Obol, is higher than 50% for 2 out of 9 AGNs. 
Correcting for incompleteness as explained in Section 6.4.1, 
we end up with the following conclusions: 

• Fraction of sources with only a starburst component 
detected: ~ 40%; 

• Fraction of sources with a dominant starburst component 
and a minor AGN component ~ 30%; 

• Fraction of sources dominated by an AGN: ~ 30%. 

Simply averaging the AGN and starburst components 
we find that the overall contribution of starbursts to the 
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Table 6. Relative AGN contribution to the L-band luminosity (a) and to the total bolometric luminosity, agpL: ^.nd L-band optical 
depth for the AGN component, r/^, for the complete sample (first 18 sources) and for the control sources (last 5 sources), as estimated 
with the model described in Section 6.4.2, using the equations in Appendix. The second estimate of the continuum L-band optical depth, 
T£^(t3,4), is based on the optical depth t^a of the 3.4 fim absorption feature, according to a standard extinction curve (Pendleton et 
al. 1994). Quoted errors are due to measurement uncertainties only. Systematic effects maybe relevant in ctBOL (of the same order as 
statistical ones for AGN-dominated sources, and 2-3 times larger for starburst-dominated sources). 



luminosity of ULIRGs is ~ 75%, while the AGN contribution 
is the remaining 25%. Considering that the typical errors in 
estimating a in Fig. 5 vary from a few percent to ~ 50%, and 
the propagation of errors in the computation of the average, 
we estimate the statistical uncertainty on the total averages 
to be of the order of 5-10%. 

However, the determination of obol has several caveats 
which must be taken into account. While the main uncer- 
tainties in the determination of a and r are due to mea- 
surement errors, other systematic effects may be relevant, 
or even dominant. In particular, two possible complications 
must be considered: 

(i) A completely hidden component, either AGN or star- 
burst, due to heavy obscuration in the L-band, cannot be 
excluded. This would imply a contribution to the bolometric 
luminosity from the FIR, with no indications in the L-band. 
Such a scenario has been proposed for instance by Fischer 
il'i) and Imanishi, Dudley, & Maloney 1^3) explain the 
dependence of the PAH-to-FIR flux ratio with luminosity 
in starbursts: lower ratios are observed in higher luminosity 
sources, suggesting the presence of an obscured component. 
In order to test this scenario, we compared the total lu- 
minosity predicted with our model (and therefore based on 
the observed L-band components) with the FIR luminosities 
measured by IRAS. The results, plotted in Fig. 6 show that 
in most cases the two values match, suggesting that hidden 
components do not represent an important contribution to 
the bolometric luminosity. However, several points are sig- 
nificantly off the 1:1 relation. In these cases a contribution 
from a completely obscured component cannot be ruled out. 
Moreover, while the bolometric correction used for AGN is 
independently determined (Elvis et al. 1994), the L-band to 
FIR ratio for starbursts is inferred from our sample (Fig. 3b). 
Therefore, if a completely obscured component is responsible 
of a fraction of the FIR emission, we are implicitely assum- 
ing that it consists of a pure starburst. As a consequence, 
our estimates of the AGN contribution could actually be 
lower limits in some cases. 

(ii) The spread in the L-band to bolometric ratio, both 
in starbursts and in AGNs (Fig. 3b) makes the determina- 



tion of the parameter K rather uncertain. Moreover, the 
assumed average ratio for AGNs is that obtained from a 
sample of quasars (Elvis et al. ITll) which could be signifi- 
cantly different from that of AGN-dominated ULIRGs, due 
to possible different amounts of circumnuclear dust. Taking 
into account this uncertainty, we conservatively estimate a 
range for K from ~ 30 to ~ 300. By propagating the error 
on Equation 5 we obtain that this is the dominant uncer- 
tainty on aBOL, at least for the sources not dominated by 
the AGN contribution in the L-band. For example, in the 
cases of NGC 6240 and IRAS 23128-5919S (Tab. 6) we ob- 
tain aBOL(K = 30)=3%, UBOhiK = 300)=0.3%. Note how- 
ever that even assuming such large uncertainties, the basic 
result that in these sources the AGN total contribution is 
negligible still holds. We note that our result is not in agree- 
ment with that of Lutz et al. who estimate an AGN 
contribution in the range 25-50%. This complex source will 
be discussed in more detail in a forthcoming paper (Risaliti 
et aLB. 

Summarizing the above analysis, we conclude that our 
simple model provides a reliable estimate of the relative 
AGN/starburst contribution for statistical purposes, with 
possible exceptions in single sources due to heavily obscured 
AGN or starburst components. 

6.3.^ Comparison with previous results 

Restricting our analysis on the complete sample made of the 
ULIRGs listed in Table 6 (14 sources, 4 with resolved double 
nuclei) we note that an AGN was found in 9 cases. In two 
cases the AGN is the dominant contribution to the bolo- 
metric luminosity, while in the other 7 cases the starburst 
is dominant. In the following we compare these results with 
the classification provided by studies at other wavelengths 
(Table 5): 

• X-rays: the classification is in all cases in agreement 
with the results of X-ray spectroscopy. 

• Mid-IR: the mid-IR analysis of G98, based on coro- 
nal lines and PAH emission features in ISO spectra, is 
in agreement with our conclusions for the five starburst 
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Figure 6. Comparison between the total infrared luminosity in- 
ferred with our model from the observed L-band emission, assum- 
ing standard bolometric corrections (see text for details), and the 
luminosity measured by IRAS. 



ULIRGs. Signatures of an AGN are clearly found in two 
sources (MKN 231, MKN 273). One source (IRAS 05189- 
2524) is not in the G98 sample. One source (IRAS 19254- 
7245) show marginal evidence of the presence of an AGN (its 
position in the diagnostic diagrams is "half way" the star- 
burst and AGN regions, but a reanalysis of Charmandaris et 
al. (2003) show clear evidence of the AGN). The remaining 
four sources (IRAS 20551-4250, UGC 5101, IRAS 23128- 
5919, NGC 6240) are classified as pure starbursts. 

• Optical/Near-IR: spectral evidence of the presence 
of an AGN was found in four sources. All other objects have 
a starburst/LINER classification. Among them, the objects 
with an L-band detected AGN are NGC 6240, IRAS 20551- 
4250, IRAS 23128-5919, UGC 5101). 

Summarizing, the above comparison shows that only L-band 
and hard X-ray spectroscopy is effective in finding evidence 
of AGNs in starburst-dominated and/or heavily absorbed 
ULIRGs, while the analysis at other wavelength fails in sev- 
eral cases. 

This result can be easily explained comparing the frac- 
tions of L-band luminosity due to AGN, a, and the fraction 
of bolometric luminosity due to AGN, ubol (Table 6): in 
ALL cases where an AGN is present, its contribtuion to the 
L-band emission is dominant, even if this is not the case for 
the bolometric emission in 7 out of 9 cases. The same con- 
sideration applies to the hard X-ray band (E>2 keV), due to 
the much higher fraction of bolometric luminosity emitted 
in the X-rays by AGNs (several percent) than by starbursts 
(< 10""*), but not to any other spectral region: 
- at wavelengths between the hard X-rays and 2-5 fim the 
AGN emission is in almost all cases completely absorbed by 



dust; 

- at wavelengths higher than a few micron the AGN emission 
is diluted by the dominant starburst emission. 

Concluding, the L-band and hard X-ray spectroscopy 
are the best spectral regions in order to investigate the en- 
ergy source in ULIRGs. As pointed out in Section 6.1, the 
X-ray band analysis is limited to a handful sources with the 
current instrumentation. Therefore, L-band spectroscopy is 
the most promising way to extend our understanding of the 
AGN/starburst nature of ULIRGs in larger samples. We 
stress again that the diagnostics developed here has been 
suggested by the analysis of our high quality ISAAC spectra, 
but once calibrated does not require particularly high signal- 
to-noise data, and is capable to solve cases which were not 
easily understandable in previous analyses. A particularly 
striking example is the luminous infrared galaxy NGC 4418 
{Lin ~ 10^^'^ Lq), which show clear signatures of the pres- 
ence of an AGN in the mid-IR (Dudley & Wynn- Williams 
licl Spoon et al. [sih. but not in the L-band, according to 
Imanishi et al. 1122^ . However, when our diagnostic diagrams 
(Fig. 3a and 3b) are applied to its L-band spectrum, we find 
that it clearly lays in the AGN region. Our model indicates 
an interpretation similar to that for IRAS 20551-4250, with 
a heavily buried AGN (r ~ 3 — 4) and a starburst component 
responsible for ^^70% of the bolometric luminosity. 

6.4 Composition of the circumnuclear medium in 
ULIRGs 

We already noticed that the correction used to estimate 
the intrinsic L-band luminosity in heavily obscured AGNs, 
based on the optical depth of the 3.4 fim absorption feature, 
leads to values log 7?£°'"'" > for three out of the 4 sources 
with r3,4 > 0.3. In particular, we obtain log R"}^^^ = 2.7 for 
IRAS 19254-7245, logi?^""' = 3.2 for IRAS 08572+3915, 
logi?£°''' = 2.3 for UGC 5101. These unphysically high 
values can be discussed considering their relation (a) with 
the L-band continuum optical depth tl estimated with the 
model presented in the previous Section, and (b) with the 
gas absorption measured in the X-rays. 

(a) The high optical depth of the 3.4 ^m absorption 
features in three objects implies an L-band extinction of 
several magnitudes, assuming the standard absorption law 
of Pendleton et al. 136l) {Al ~ 12±4r3.4). Even considering 
the lowest allowed value, i.e. Al = 8x3.4, we would still 
have an L-band to total infrared ratio log > for two 

sources (IRAS 19254-7245 and IRAS 08572-^3915). 

On the other hand, the extinction tl obtained from the 
spectral decomposition presented in the previous Section is 
significantly lower in the three objects with the highest r3.4, 
and compatible with an intrinsic L-band to bolometric ratio 
log R^L^^ ~ 0. The comparison between the two estimates of 
the optical depth is shown in Table 6. 

The above results imply that the abundance of hydro- 
carbon dust grains, responsible for the 3.4 /xm absorption, 
is on average higher in these ULIRGs than in the Galactic 
ISM. This is the case for all the three sources with strong 
absorption at 3.4 ^ra {t^a ~ 0.8 — 0.9). We note that this is 
not necessarily a general feature in ULIRGs, since selection 
effects could play a role: such high r3.4 can be observed only 
in objects where the dust grains responsible of the 3.4 /xm 
absorption are overabundant, because otherwise the contin- 
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uum extinction would be so high to prevent the AGN emis- 
sion from being observable. 

(b) The column density measured in the X-rays are 
Nh > lO^"* cm-2 for UGC 5101 and IRAS 19254- 7245jbased 
on the XMM-Newton observations (Imanishi et al. |20| and 
Braito et al.|3) and Nh > W^^ cm'^ for IRAS 08572-^3915 
from a BeppoSAX observation (Risaliti et al. 2005, in prep.). 
Assuming a standard Galactic extinction curve and dust-to- 
gas ratio (Savage & Mathis 1481 we have Al ~ 1.2 Nh,23, 
where Nh,23 is the column density in units of 10'^'^ cm"'^. 
This would imply Al > 12 for UGC 5101 and IRAS 19254- 
7245, and > 120 for IRAS 08572-^3915. Again, these 
values are not acceptable, because they would imply a too 
high intrinsic continuum. We conclude that the dust-to-gas 
ratio in heavily obscured ULIRGs is at least a factor of a 
few smaller than in the Galactic ISM. This is a well known 
property of the absorbing medium in Seyfert Galaxies with 
moderate absorption, i.e. Nh < W^^ cm~'^ (Maccacaro, Per- 
ola, & Elvis 13^. Maiolino et al. l3lf) . Here we find that this is 
also the case for higher luminosity and higher column den- 
sity AGNs. 



7 CONCLUSIONS 

In this paper two works are presented: (1) L-band spectra 
obtained with ISAAC at VLT are analyzed for a sample of 7 
Ultraluminous Infrared Galaxies, and (2) a detailed discus- 
sion on L-band diagnostics for ULIRGs is performed, using 
a complete sample of local bright ULIRGs, obtained merg- 
ing our observations with the available L-band spectra of 
sources in the northern hemisphere. 

1. The spectral analysis of the sample of 7 bright south- 
ern ULIRGs show the power of ISAAC in revealing the 
nature of these objects. Six sources are known to have 
double nuclei. In four cases we obtained a spectrum for 
both nuclei. A clear signature of the presence of an AGN 
was found in three nuclei (IRAS 19254-7245S, IRAS 20551- 
4250, IRAS 23128-3915S). In two of these three objects 
(IRAS 19254-7245A, IRAS 20551-4250) an absorption fea- 
ture at 3.4 ^m, due to hydrocarbon dust absorption was de- 
tected. These two sources have also the steepest continuum 
(r > 2). Water ice absorption at 3.1 fim was also detected 
in several objects. We confirm that 3.4 absorption is a 
typical signature of an AGN, while 3.1 absorption with 
moderate optical depths (ra.i < 0.5) is found in sources with 
no indication of an active nucleus. 

The high quality of the spectra obtained with ISAAC 
demonstrates that L-band spectroscopy can be a powerful 
tool to investigate the nature of relatively faint ULIRGs. For 
example, most of the 118 sources of the 1 Jy sample (Kim 
& Sanders )24j) are probably bright enough for such a study. 
This is not true for studies at other wavelengths, such as in 
the X-rays, which are in principle extremely effective in dis- 
entangling AGNs and starbursts but can only be performed 
with the 10/15 brightest sources with the currently available 
instrumentation. 

2. We performed a complete analysis of the L-band spec- 
tral diagnostics on a sample consisting of our seven southern 
ULIRGs, 7 northern ULIRGs selected with the same crite- 
rion (60 fj,m fiux density 5*60 > 5.2 Jy) and 5 more sources 
which are known to be either AGN-dominated or starburst 



dominated, with a slightly lower bolometric luminosity. All 
these sources have a reliable, independent starburst/ AGN 
classification based on observations at other wavelengths. 
The main results are the following: 

• In all known AGN, at least one L-band AGN indicator 
is detected, while no AGN indicators are found in any object 
classified as starburst. However, no single indicator provides 
a 100% correct classification. 

• When more indicators are combined, a 100% correct 
starburst/ AGN classification is possible. We found that 
AGNs and starbursts are completely separated in two differ- 
ent bi-dimensional diagrams, the first obtained plotting the 
3.5-4 fim continuum slope, F, versus the 3.3 ^m PAH emis- 
sion feature equivalent width, EW3.3; the second plotting 
-EW3.3 versus the extinction-corrected L-band to bolometric 
emission ratio. 

• ULIRGs hosting AGNs can be clearly separated into 
two groups. The first consists of objects where the AGN is 
not heavily absorbed in the L-band. The continuum is dom- 
inated by the re-emission by hot dust of the AGN direct 
emission; the continuum slope (in a A — /a plane) is flat; the 
PAH emission feature is weak or absent; absorption features 
are absent. The second group is made by objects heavily 
absorbed in the L-band. In these cases the emission is dom- 
inated by by the reprocessing by cold dust. The continuum 
is steep (F ~ 2) due to reddening, and strong absorption 
features are present. 

• We use a simple model for the spectral deconvolution of 
the AGN and starburst components. We show that it is pos- 
sible to estimate the relative AGN and starburst contribu- 
tions to the bolometric luminosity through L-band spectral 
indicators. 

• The application of our model to a representative sam- 
ple of ULIRGs in the local Universe shows that AGNs are 
present in ~ 60% of ULIRGs. In about half of these cases the 
AGN contribution to the bolometric luminosity is dominant. 
Overall, the AGN contribution to the ULIRG luminosity in 
the local Universe is ~ 30%, the remaining 70% being due 
to starburst emission. 

• The comparison between gas absorption (from X-ray 
observations), continuum extinction in the L-band, and the 
narrow absorption feature at 3.4 jj,m suggest that heavily ob- 
scured ULIRGs have a low dust-to-gas ratio and a different 
extinction curve with respect to the Galactic ISM. In partic- 
ular, hydrocarbon dust grains are overabundant in three out 
of four cases with respect to the Galactic dust composition. 
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APPENDIX A: ANALYSIS OF THE AGN/SB 
DECONVOLUTION MODEL 

The relations plotted in Fig. 5, and expressed in Equations 3 
and 4 can be inverted in order to work out the explicit de- 
pendence of the two model parameters, a and tl, on the 
observed parameters, EW and F: 

EWsbT - EW 



0{EWsB - EW' 

EWsB - EW 



(Al) 



(A2) 



{EWsB - EW) + EW * e-'L 

where /3 = 1.75 is the slope of the extinction curve in 
the L-band, EWsB is the intrinsic equivalent width of the 
3.3 /im PAH feature in a pure starburst spectrum. 

From these two equations we obtain the errors on tl 
and a: 
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where A[EW) and A(r) are the errors on our measure- 
ments. 

These equations have been used to estimate the errors 
in Table 6. Moreover, they can be analyzed in order to quan- 
titatively estimate in which regions of the (EW, V) parame- 
ter space we can measure the relative SB/AGN contribution 
with good precision. In Fig. Al we plot the relative errors, 
A{tl)/tl and A{EW)/EW versus the measured EW , for 
different values of F, and assuming A{tl)/tl = 0.1 and 
A{EW)/EW=0.2. 

The regions with a bad determination of a (with errors 
higher than 50%) are those with high equivalent width of the 
3.3 jim emission feature and a steep spectrum, correspond- 
ing to the upper right region in Fig. 5 where the model lines 
get closer while approaching the pure starburst point. This 
is also clear from Eq. A4, where the value of Aa diverges 
for EW EWsB. 
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